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The linear approximation of scalar-tensor theories of gravity is obtained in the physical (Jordan) 
frame under the 4+0 (covariant) and 3+1 formalisms. Then the weak-field limit is analyzed and the 
conditions leading to significant deviations of the 1/r 2 Newton's law of gravitation are discussed. 
Finally, the scalarization effects induced by these theories in extended objects are confronted within 
the weak-field limit. 
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I. INTRODUCTION 



In the last decade there has been an increasing interest 
in the so-called scalar-tensor theories of gravity (STT; see 
Rcf. for a review) in view of the possible deviations 
that Einstein's general theory of relativity (GR) could 
show in the framework of several upcoming observations. 
Among these, there has been a considerable interest in 
detecting a scalar-wave component (spin-0 waves) in ad- 
dition to the ordinary gravitational waves (spin- 2 waves) 
predicted by GR (see for instance (|, ||, ||, ||, |(| ) . Among 
the potential emitters of such scalar-waves are the com- 
pact binary systems and neutron stars. Yet, the binary 
pulsar has showed no appreciable deviations from GR due 
to the emission of scalar waves, and then one is only able 
to put stringent bound in the STT parameters or cou- 
plings JtJ. Moreover, through the effects of spontaneous 
scalarization in neutron stars (?], [|, fLO] , it seems that 
scalar-waves have also a small chance to be observed in 
detectors like VIRGO or LIGO within a few hundreds of 
kiloparsecs ||. However, there is the hope that resonant 
mass detectors of spherical shape or interferometers like 
LISA can resolve the existence of a scalar-gravitational 
wave |,!|. 

At the large scale (cosmological scales), STT have been 
proposed as models for dark energy that can replace the 
cosmological constant fTlJ| , and also for explaining some 
other features of the galaxy distribution in our universe 
fL2| . However, detailed analysis show that when all the 
observational constraints (both cosmological and local) 
are taken into account, the simpler STT fail or are ex- 
tremely constrained At this point one cannot only 
but recognizes the predictive power of GR, confirmed by 
experiments and observations that were maybe far from 
being imagined by Einstein at the time of the creation of 



his theory. At smaller scales (scales of the order of me- 
ters or kilometers) the spurious discovery of a fifth force 
renewed the idea about the existence of new fundamental 
fields of meter-range |l4|]. At this regard, considerable ef- 
fort has been put in measuring gravitationally such kind 
of interactions that can mimic a varying gravitational 
"constant" G. 

In this report the linear limit of STT is reviewed under 
the usual 4+0 formalism. The weak-field limit is taken 
and confronted in view of the past and recent experi- 
ments (ground based and satellite mission) intended to 
test the gravitational interaction between bodies. The 
analysis is performed in the Jordan frame, and there- 
fore no intermediate unphysical variables or transforma- 
tions are introduced. A Yukawa potential appears as a 
new term in addition to the ordinary Newtonian poten- 
tial due to the non-minimal coupling between the scalar 
field with gravity, the entire gravitational potential hav- 
ing been identified with the potential of test particles. It 
is then argued that the Yukawa potential can produce 
significative deviations on the 1 /r 2 Newton's gravitation 
law if the scalar "particles" are massive enough. 

An appendix analyzing the linear limit of STT under 
the 3+1 formulation of GR is to be found at the end. 



II. SCALAR-TENSOR THEORIES OF GRAVITY 

The general action for a scalar-tensor theory of gravity 
with a single scalar field is given by 



x^/^gd 4 x + S\gp V ,ip] 
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where ip represents collectively the matter fields (fields 
other than <fr). 

The representation of the scalar-tensor theories given 
by Eq. (M is called the Jordan frame representation. One 
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can parametrize the same theories as 
1 



where 
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§R - ^^(V$) 2 + 2$A($)| 



(2) 



wbd(*) 



A($) := - 



m > 

8ttG $ 



(3) 
(4) 

(5) 



For instance, the Jordan-frame representation of the 
Brans-Dicke theory with to =const. corresponds to 
F = 2nGo(j) 2 /uj and V((f>) = 0. It is also customary 
to parametrize the scalar-tensor theories in the so-called 
Einstein frame by introducing non-physical fields as fol- 
lows, 
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F*(</>*) = F(4>) . 
so that the action Eq. (|l|) takes the form 
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, (?) 
(8) 



S[g*4>*A] = 
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[R* - 2(V*0*) 2 - V*{(j}*)} 



Xy^d 4 x + S[g* /F* («/>*) 



(9) 



where all quantities with '*' are computed with the non- 
physical metric and 0*. 



We emphasize that although the equations of motions 
obtained from the Einstein frame are simpler than those 
from the Jordan frame, in the sense that the field (f>* ap- 
pears to be coupled minimally to the non-physical metric, 
the matter equations derived from the Bianchi identi- 
ties V*G* V = will have sources, i.e., V^T* ^ ^ 0, 



where here T* " u = T^/F* 



) is the non physical 
energy-momentum tensor of the matter fields ip. How- 
ever, in the Jordan frame the matter equations resulting 
from the Bianchi identities V M G Miy = turn to satisfy 
V^T^" = 0, reflecting explicitly the fulfillment of the 
Einstein's weak equivalence principle (that is the origin 
of the name "physical metric" ) . 

In the following the Jordan frame representation of the 
scalar-tensor theories will only be used. The equations 
of motion obtained from a variational principle using Eq. 
©are 
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(ii) 
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(12) 
, (13) 
(14) 



Fd*V - 2{d 4) F)V - \(d+F) (l + S?) (V0) 2 + ±(d F)T matt 



D ^ = jjMEZ ■ (15) 



where T ma tt stands for the trace of T^ att and the sub- 
script "matt" refers to the matter fields other that <\>. 
Now, the Bianchi identities imply 

V M T^ = . (16) 

However, the use of the equations of motion leads as men- 
tioned to the energy-conservation equations of matter 

V„T^ tt = , (17) 

which implies in the case of test particles, that bodies 
are subject to no other long range forces than the gravi- 



tational ones (free falling particles). In other words, the 
scalar field being not directly coupled to matter no direct 
interaction between the scalar field and matter arises. 
The scalar field (f> will only interact with matter gravita- 
tionally, i.e., only through the curvature effects. 

The final form of the field equations have exactly 
the same form as in general relativity with an effective 
energy-momentum tensor. This means that one can for- 
mulate the Cauchy problem for scalar-tensor theories ex- 
actly in the same manner as in GR (see the Appendix). 
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III. LINEAR LIMIT OF STT 

The linear limit of STT has been analyzed in the past 
by many authors. Wagoner Q was one of the first in an- 
alyzing the gravitational wave emission in STT and the 
weak-field approximation (an updated analysis was per- 
formed in Q). Such an analysis was performed in the 
Einstein frame. More recently, Pimentel and Obregon 
fl5| performed a similar analysis in the Brans-Dicke rep- 
resentation as in Eq. (||). 

As mentioned, the linear limit of STT treated here is 
performed in the Jordan frame. As usual, we consider 
first order perturbations of the Minkowski spacetime: 
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where e< 1 and the knott indicates quantities at zero or- 
der. In the 4+0 covariant formulation one can introduce 
the combination 



1 



(26) 
(27) 



where 7 = 7 M M and n is a gauge constant to be fixed 
later in order to simplify the equations. The resulting 
linearized Einstein equations are 



= 8%G o f^ + K(dlJ-ri„vO n 0) . (28) 

where T° = results from the self-consistency of 
the perturbations at first order. Here is the 
D'Alambertian operator compatible with the flat met- 
ric rjfiu- Moreover, the flat metric is used to raise and 
lower indices of first order tensorial quantities. 
The Lorentz gauge 
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(29) 



The linear approximation of the effective energy- 
momentum tensor Eq. ( jll| ) and the Klein-Gordon Eq. 
(Ill) turn to be 
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where we used the following conditions 

v(<t> ) = = vtf = , 



(31) 
(32) 
(33) 

(34) 
(35) 



resulting from the consistency at first order of the lin- 
earized Einstein and Klein-Gordon equations and assum- 
ing that T^ att = 0. In this way the wave equation Eq. 
(|30|) becomes 
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FoJ \ ^ 



4> - r)v,v a r)(/> 
(36) 



Note that at this order the application of the ordinary 
divergence in Eq. (^) and the use of the Lorentz gauge 
Eq. ( |29| ) leads to the energy-conservation of the matter 
perturbations: d^T^f" = 0. 
The choice 



(37) 



for the gauge constant simplifies the Eq. ([36|) consider- 
ably. So summarizing, we have the following wave equa- 
tions for the gravitational and scalar modes: 



□ 



16nG ~ matt 

Fn 



(38) 
(39) 
(40) 



with the constants Wq, a and k given by Eqs. (|3|), ( |34| ) 
and (|37]) respecteively. The analysis of propagation of 
gravitational and scalar waves will be not pursued here, 
this has been done elsewhere (see Refs. [0, g[ 0, 0, § ). 



generalizes the usual Lorentz gauge of GR, and can be 
imposed to simplify the equations. Then, from Eq. (I2F 
the resulting wave equation is 



□ 



(30) 



A. The weak-field approximation 

In the weak-field approximation one considers slow 
varying fields and sources such that the temporal deriva- 
tives in the D'Alambertian can be neglected and T'™ att w 
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pS^S®. The wave equations read then, 



3A- _ 167rG 0y, mat t 



>A% V = — 



Fn 



3 A0 - ml4> = 47ra^f matt 



(41) 
(42) 



where 3 A is the three dimensional Euclidean Laplacian. 



The condition for slow varying sources T^ &tt 
T™ att , leads then to the following equations 







3 ~ 16ttGo~ 


(43) 


pi 

3 A</> - ml4> = -4na-£p , 


(44) 


7oi = const. = , 


(45) 



where p := T^ u = -T matt . The constant solutions for 
7 Mi arise by demanding a well behavior of the metric at 
spatial infinity (without lost of generality the constants 
can be gauged out). Thus, with jij = and Eqs. (|26| ) 
and @, it yields 



%v = 7/^ - Vij.i< (loo + 2/c<^ 



(46) 



and therefore 



(700 + 



700 = 2 700 + K 



(47) 



Since the geodesic equation for slow particles leads to the 
identification of 



N 



2 700 



with the Newtonian potential p6[ , it turns 

1- k 7 



■4700+2' 



(48) 



(49) 



The two gravitational degrees of freedom that remain in 
the weak-field limit are thus 700 and 4> (joj = 0, and 
jij is given in terms of 700 and 4>)- 

The solutions of Eqs. (43) and (0) are respectively 
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Fq j \x — x' 
F^ f p(x')e- ma \ s - s '\ 
( FoJ \x-x'\ 



where -0 : = — t7oo- So finally 
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(50) 
B.C. , (51) 
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\x — x' 



hi ^'\p- m o\S-x 



B.C. 



(52) 



It is important to mention that the details of the general 
solution (interior plus exterior solutions, both matched 
continuously at the surface of the extended object) de- 
pend strongly on the boundary conditions. 



and asymptotic flatness conditions 



lim r 
lim r 



, 
. 



(54) 
(55) 



B. Regular solutions 

Explicit solutions of the gravitational and scalar fields 
depend on p(x), and as mentioned on the boundary con- 
ditions as well. Clearly, a detailed solution of the whole 
system involves also the equation of hydrostatic equilib- 
rium and the equation of state of matter Jl7f| . However, 
here we do not want to be so explicit and only exhibit 
the possible deviations that the Newton's gravitation law 
can suffer depending on the values of the parameters. For 
simplicity spherical symmetry is assumed. Wc seek for 
solutions for tp and <j) with regularity conditions at r = 



d r ip = = 9 r 



(53) 



The solutions arc then given by 



ip(r) 



+oc m(f) 
F b J r ~f 



df , 



where 



(56) 
(57) 



m(r) := 4ir I p(f)f df , 
Jo 

q(r) := 4tt / ( a^-p(f) - -^0 I f' 2 df 

J \ F Q 4tt 



(58) 
(59) 
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One can define global quantities as 



M 



= 4-7T 



p{f)f 2 df 
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a Y p{r) ~ iT 



(60) 



(61) 



Note that in the minimal-coupling case (Fq = 0), Q = 0, 
since the only regular solution of the Helmholtz equation 
is <f> = 0). 

If the density p has compact support (as it is usually 
the case for astrophysical bodies), then it vanishes for 
r > R, where R is the radius of the body. The solutions 
for ip and <f> matched continuously at r = R can be then 
written as 



ip(r) = 



■^-^[ fi ^rfr for r<R 

RFq Fo Jr r 2 — 



Ggh 
F r 

Q 



for r > R 



S(ray + / r i? #df for r<R 



r{l+m R) 



where 



m(r) 
q(r) 



r r<R 

4-7T / p(f)f 2 df 





(62) 
(63) 

(64) 



<j) f 2 df , (65) 



and M := m(R), Q := 



IV. OBSERVATIONAL CONSTRAINTS 

In the following we consider only the exterior solutions 
and confronted them with observations. Therefore, the 
exterior Newtonian potential is given by 
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2Gq -Fq 



(66) 



where cr is a constant that depends on the global prop- 
erties of the source. Namely, 

Qe moR F 



(l + m R)MaF^ 

p ra R ( m 2p o 



1 + moR 
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MaFa 



(r)r 2 dr 



(67) 



It is customary to express the coefficients involving Fq in 
terms of the effective Brans-Diecke parameter Eq.([|) 



InGpFo 



Then, 
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4tt (3 



2w° D ) 



2o;^ d Vq' 
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and therefore 



G M 
F r 



3 + 2^ D 



(68) 

(69) 
(70) 

(71) 



Recently experimental bounds on the strength and 
range of a Yukawa potential that could arise from "fifth 
force fields" have been analyzed in two kinds of experi- 
ments. The first kind consists in the analysis of gravi- 
tational signals induced by variations on the mass of a 
lake |18|| . Such experiments probe basically fields with 
a range A from meters to some kilometers and strength 
\j3\ € [10 _4 ,10 2 ]. Here the coefficients corresponds to 
P = ae- r / x /(3 + 2^L) and A = l/m (cf. the curve 
±/3 = ±/?(A) in Ref.[n8[) . The second kind of experi- 
ments probes variations of the Newton's gravitation law 
at scales of two Earth's radii by measuring the gravita- 
tional effects on the orbit of the laser-ranged LAGEOS 
satellite Assuming that r <C A, so that effects of 

an "intermediate-range" force are taken only to order 
(r/A) 2 , it turns that |/3| < 10 -5 - 10~ 8 . These bounds 
are even more restrictive than those from the Earth based 
experiments quoted above. 

At solar-system scales, Viking-like experiments restrict 
w bd > ^cxp, where cj C xp ~ 3000 corresponds to the cur- 
rent lower bound on Wg D p0| . This bound results from 
the parametrized post-Newtonian(PPN) approximation 
and translates into a bound on (3 and A: 



< 
A < 



-r/A 



(3 + 2w cxp ) 

/ 3 + 2oj oxp 
V 2v exp V " 



1.66 10~ 4 fi e- r/x . (72) 
' (73) 



Jo 



If we assume for instance that V — m 2 cf) 2 and mo <C eV 
, for example, too ~ 10~ 21 eV, then it turns that A < 13 
kpc, and moreover from Eq. ( |67| ) one can expect that 
cr ~ 1. Then, the most stringent upper bound ~ 1/(3 + 
2djg D ), imposed on /3by the solar system experiments 
arises for long range scalar fields A ~ kpc (i.e., galactic- 
scale ranges) : 
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(3 < 1.66 10~ 4 
r <C kpc . 



(74) 
(75) 



However, at the weak-field limit, a violation of the con- 
dition ( |74| ) do not imply significative violations of the 
1/r 2 gravitational force but rather a renormalization of 
the gravitational constant or of the self-gravitating mass 
(see below). At the PPN level, the effects of such "mass- 
less" scalar fields manifest sensitively in the propagation 
of electromagnetic radiation (radar-echo delay) . 

Conversely, the solar system experiments impose weak 
bounds on scalar-tensor theories if the scalar field has 
short range, since then the Yukawa interaction would 
not be felt out of the self-gravitating body. However, as 
pointed out above, at regimes of Earth range the "fifth 
force tests" start playing some role and therefore devia- 
tions of the 1/r 2 gravitational force could arise. More- 
over, if A ~ kpc, then also at galactic scales a violation 
of the 1/r 2 law could be expected. 

Effective gravitational constant and total gravi- 
tational mass. If we assume for instance that V — m 2 </> 2 
then the self-consistency at first order according to the 
conditions ([35|), implies that m = or 4>o = (or both). 
In order to facilitate the interpretation of the following 
considerations, let us assume that V = which corre- 
sponds to a massless scalar field. 

Now, the usual interpretation of fundamental con- 
stants in physics is by identifying the coefficients of some 
fundamental fields in the corresponding Lagrangian. For 
instance, one identifies the electric charge from the coef- 
ficient of 1/ ' e 2 F llv F llv in the electromagnetic Lagrangian, 
or the mass of fermions as the coefficient of mifji/j m the 
Dirac Lagrangian. In the same way, the effective gravita- 
tional constant can be identified as the coefficient of the 
Ricci scalar in the Lagrangian (0) , that is by Eq. (14). 
This means that at zero order [Ell, 



^eff — TT 
-TO 



(76) 



So rather to define an effective gravitational constant 
fromEq. @: 



G := 



2G 
F 



3 + 2^ D 



(77) 



as it is usual in most treatments that have analyzed the 
weak field limit, we shall keep the definition of Eq. (76). 
Instead one can define an effective mass 



M 



oil' 



M 



= M 



1 



2G 

HQ 
2G ' 



2M 



(78) 



At this regard one mentions that it is the combination 
GoffMeff that always appear in the Newtonian potential 
( [FT| ) and not the separate quantities. So whether one 
adopts one or the other definition for G e ff or M e g, the 
combination of both provides the same number. 

Naively one could then expect that M e ff is the total 
gravitational mass of the system. However, this is not 
the case, since the total ADM -mass formula in the weak- 
field limit gives [see Eq.( A64 ) below] 



DM 



f 00 d 3 x 



(79) 



which turns to be 



-Madm 



M 



= M — 



(Hf 
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= 2M 



2w BD 



HQ 

2G 



(80) 



This result is not surprising if one takes into account that 
it is in fact the active mass M c ff and not the ADM-mass 
which is measured by orbiting test particles B9L 

It is customary in the literature of gravitational physics 
to introduce different mass definitions of a body. For in- 
stance, tensor, scalar and inertial (or sometimes referred 
to also as Keplerian) masses Mr, Ms and Mj 
In the weak-field limit these arc given by 



M T 
M s 
Mi 



M 



M 



3 

M eff 



2w BD 



HQ 

2G 



(81) 
(82) 
(83) 



In this limit, the tensor mass corresponds to the total rest 
mass since M is just the integral of the energy-density 
of the self-gravitating body which at first order coincides 
with the integral of the rest-mass density over the proper- 
volume (2^] . The scalar mass ( |82|) is simply proportional 
to the scalar charge of the body. The ADM-mass then 
reads in this limit as, 



Madm = M T - M s = Mi - 2M S = 2M T - Mi 



(84) 



and the other masses can be expressed in terms of the 
following combinations: 



Mi 



Madm 
Madm 



M s = M A 



M s = Mi - M s 
2M S = M T + M s 
1 



(85) 
(86) 



DM 



2M^ 



s 



(Mr -Madm) • (87) 



Spontaneous scalarization. The phenomenon of 
spontaneous scalarization J?], [| ||, arises in scalar- 

tensor theories within a dense compact object (e.g. neu- 
tron star). This phenomenon corresponds to a "sud- 
den" appearance of a non-trivial configuration of the 



7 



scalar field when the extended object turns to be com- 
pact enough and in the absence of scalar-field sources (as 
opposed to induced scalarization) . In other words, the 
spontaneous scalarization phenomenon in neutron stars 
corresponds to a solution of field equations with fixed to- 
tal baryon number and a zero value for the asymptotic 
boundary condition 4>q. In this way there are two possi- 
ble configurations: one with a trivial scalar field (0 = 0) 
with larger ADM-mass and another with a non-trivial 
scalar field with lower ADM-mass, both configurations 
having the same total baryon number and 0o = 0. It 
has been argued that spontaneous scalarization is a non- 
perturbative effect || and therefore it is not expected 
to appear in the weak-field limit. This conclusion can be 
checked by analyzing the regular solution for (f> in Eq.(44). 
In order to have spontaneous scalarization one would re- 
quire Fq = and at the same time Fq ^ (with </>o = 0) , 
both conditions needed to obtain a trivial and a regu- 
lar non-trivial configurations for <fi. Since a multi-valued 
derivative is impossible for analytic functions F(<f>), one 
concludes that spontaneous scalarization takes no place 
at the weak-field level. In other words, in the weak-field 
limit, the asymptotic condition O = fixes for one and 
for all the value of the constant Fq . This value has to be 
zero for the (f> = configuration to exist (the ADM-mass 
( [sol ) has then the unique possible value Madm = M) 
and different from zero to find the non-trivial scalariza- 
tion solution <fr ^ (both conditions are then impossi- 
ble). However, in the non-perturbative analysis and still 
in the massless case V = 0, one concludes from Eq. jl5| ) 
that a necessary condition for spontaneous scalarization 
to ensue is that F has a critical point in <fi = (i.e., 
d$F\$ = o) — 0, with Tmatt ^ 0), since only in that case 
the solution <f) — is possible. For instance, an even func- 
tion F(<j>) verifies that condition. The other necessary 
and sufficient condition to obtain spontaneous scalariza- 
tion in the full theory results by proving that there exists 
a non-trivial stationary regular solution of the non-linear 
Klein-Gordon equation Eq. ( |i~5| ) for a configuration leav- 
ing the rest-mass unchanged and for the same asymptotic 
condition 4>q — 0. This has been done only numerically 
and for some class of STT ||, ^j, For example, the 
particular case F = 1 + 167r£</> 2 , was analyzed non per- 
turbatively in Ref. JlO| via a numerical analysis and the 
conclusion was that only within highly compact neutrons 
stars (strong field regime) the phenomenon appeared. 

As opposed to spontaneous-scalarization there is the 
induced scalarization. This arises when <f>Q ^ 0, and since 
any analytic (non-trivial) function F(4>) implies Fq =/= 
for cj)Q ^ 0, it has a weak-field limit analogue (it would 
be rather artificial to expect a Lagrangian which depends 
explicitly on 4>q and moreover that Fq = 0). The simplest 
non-trivial example is F = l+Wir^cf). In this case F'L = 
const, and therefore a trivial solution for <f) is not allowed 
[cf. Eq.© or ( ||) with V = 0]. An artificial coupling 



F = 1 + 167r£(0 - </> ) 2 leads however to F'\<f, = 0. In 
other words, with the asymptotic condition cf>Q ^ and 
with a coupling function F{<f>,4>o) such that d$ F = 0, 
it is impossible to have a trivial scalar field in any limit 
(weak or strong regime), since a trivial scalar-field <f> — 
const, would not be a solution of Eq. ( |l5| ) for Tmatt 7^ 0. 
Induced scalarization then always ensues. The fact that 
the induced scalarization is present in the weak-field limit 
produces a smooth transition to the strong field regime 
(as opposed to spontaneous scalarization). Namely, the 
curve Qus.Mbar (scalar charge vs. total baryon mass) 
do not show any discontinuity contrary to spontaneous 
scalarization (cf . Ref. Jj], ||, [l(J ) . 



V. CONCLUSIONS 

The linear limit of scalar-tensor theories of gravity 
leads to the prediction of a spin-zero gravitational mode 
(scalar waves) in addition to the well known spin-2 
modes. In the weak-field limit, STT predicts an effec- 
tive gravitational potential with a Yukawa contribution 
for massive scalar fields. However, for ultra-light scalar 
fields (scalar fields of galactic or cosmological range) the 
Yukawa term is highly supressed as far as solar sys- 
tem experiments are concerned. Thus, in such a case 
the effective gravitational potential leads essentially to a 
Newtonian gravitational force ~ 1/r 2 with an effective 
mass coefficient which is proportional to the rest-mass 
of the self-gravitating body and which depends explic- 
ity on the effective Brans-Dicke parameter. This latter is 
highly constrained by the post-Newtonian approximation 
via Viking-like experiments. Therefore, even if a STT 
does not produce a deviation of the 1/r 2 Newton's law, 
it can however violate dramatically the post-Newtonian 
bounds. On the other hand, if the scalar field has "inter- 
mediate" ranges (meters to kilometers or even thousands 
of kilometers), then the Yukawa term leads to a "fifth 
force" field that produces deviations of the 1/r 2 New- 
ton's law and which are severely constrained by Earth 
based experiments or by satellites, notably, by varying- 
mass experiments and laser-ranged missions. STT also 
predicts the phenomenon of spontaneous scalarization in 
compact objects. For the phenomenon to take place, 
a STT requires the existence of two possible configu- 
rations, one of which corresponds to the trivial <f> = 
solution and another one that is not trivial. Both config- 
urations having the same rest-mass and the null asymp- 
totic value (j)Q = 0. The phenomenon will not be present 
in the weak-field limit since it would require that the 
non-minimal coupling function F(<f>) has a two valued- 
derivatives: i* 1 ' (0)| Q =o = and -F'(0)|<£ o =o 7^ in order 
two obtain the two possible configurations for the per- 
turbation (f>. The conclusion is that no-analytic function 
F{<f) can fullfill both conditions. As oppposed to sponta- 
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neous scalarization, induced scalarization takes place also 
in the weak- field limit, since for physical couplings one 
has F'(<f)o) 7^ (i.e. physical coupling functions F will 
not depend explicitly on <j> ), and then only non-trivial 
regular solutions are admitted. 
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This means that 

e = E i > J ^ = 2 J ^ ' S " V =J2 s r ■ (A6) 

i i i 

The projection of Einstein equations = 
4nGo (2T M „ — T a a g^ u ) in the directions tangent and or- 
thogonal to S t , followed by the use of the Gauss-Codazzi- 
Mainardi equations leads to the 3+1 form of Einstein 
equations: 

3 R + K 2 - /\,,A'' = 16nG E , (A7) 
known as the Hamiltonian constraint. 



APPENDIX A: THE 3+1 FORMULATION OF 
GENERAL RELATIVITY 



Let us consider the 3+1 or Adison-Dcser-Misner 
(ADM) formulation of general relativity in which the 
spacetime (considered to be globally hyperbolic) is foli- 
ated by a family of spacelike hypersurfaces S t . We shall 
not enter into the details of the derivation of the 3+1 
equations (see Refs. |2(| |28|, |29|). The sign con- 
vention for the 3+1 splitting of the metric is as follows: 



ds 2 = -(N 2 - N l Ni)dt 2 - 2N,dtdx l 



+ hijdx 1 dx J 



(Al) 



The extrinsic curvature of the embedings S( is given by 



-V»rij 
1 



-AT*- 

l 3 



dhi 



2N V dt 



ViiVj , (A2) 



where 3 V-,- stands for the covariant derivative compatible 
with the three-metric hij. This is to be regarded as an 
evolution equation for hij. The trace of the extrinsic 
curvature will be denoted by, 

K 



(A3) 



The orthogonal decomposition of the energy- 
momentum tensor in components tangent and orthogonal 
to £ t leads to || : 



(A4) 



where is the normal to S t . The tensor S^" is sym- 
metric and often called the tensor of constraints; is 
the momentum density vector and E is the total energy 
density measured by the observer orthogonal to S( . 

As in the 4+0 formalism, will be the total energy- 
momentum tensor of matter which can be composed by 
the contribution of different types of sources: 



(A5) 



3 V;-fsT' 4 - 3 V,A- = 8ttG Ji , 

known as the momentum constraint equations. 
Finally, the dynamic Einstein equations read 



(A8) 



d t K t j + N%K x j 



K\djN l 



K jdiN 1 



VjN- R 1 jN — NKK 1 j 



+ 3 V' 3 

= 47rG Y [{S-E)5\ - 2S 



(A9) 



3 3\ 

where S = S\ is the trace of the tensor of constraints, and 
all the quantities written with a '3' index refer to those 
computed with the three-metric hij. Moreover, under 
the 3+1 formalism tensor quantities tangent to T, t use 
the three-metri c to raise and lower their spatial indices. 
The equations (|X|) and (|AS| ) are the set of the Cauchy- 
initial-data evolution equations f or t he gra vita tional field 
subject to the constraints Eqs. (A7) and (|A8|). 

We can write an evolution equation for the trace K by 
taking the trace in Eq. (|A9|): 

3 AN -N( 3 R + K 2 ) = 4-kGqN [S - 3E] 

(A10) 



d t K + N l d t K 



where 3 A stands for the Laplacian operator compatible 
with ha. 



This can be simplified by using Eq.(A7) to give 



d t K + N l diK + 3 AN - NK tj K ij = AttG N [S + E] . 

(All) 

For our purposes it will no be necessary to write the 
3+1 equations for the matter and the scalar field. 



1. Linearized equations 

In order to linearize the 3+1 equations we assume first 
order deviations of the 3+1 metric with respect to the 
Lorentz metric n^ as follows 



N 

Ni 
hij 



1 + eN 

eN, , 
Sij + eh 



(A12) 
(A13) 
(A14) 
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At first order h lJ w 5ij — ehij in order to satisfy the in- 
verse condition huh J — 6 l j. Therefore, covariant and 
contravariant components of tensorial quantities having 
no zero order terms are identical to each other at first or- 
der. For instance N 1 = h ll Ni w eiV,-. In order to compare 
with the 4+0 (full covariant) linear approximation 



it turns out that 



700 = -2N , 
loi = -Ni , 
7ij = hij ■ 
The 3-Christoffel symbols turn to be 

1 

2 



e 3 f l 



3^2 



djijk 



dkhij + d-jhk 



(A15) 

(A16) 
(A17) 
(A18) 

(A19) 
(A20) 



Therefore 3-covariant derivatives of 3-tensors having no 
zero order terms become at first order orcfinary deriva- 



tives. For instance 3 S7 Ni ps ed 3 Ni. Then Eq.([A2j) leads 



to 



Ki 



■thij 



2d {l N £ 



with 



K' 



j 
K 

K 



1 

~2 



» eKij , 
eK , 

d t h + 2d l Ni 



(A21) 
(A22) 

(A23) 
(A24) 

(A25) 



where h := hu is the trace of the 3-metric perturba- 
tion. Using the definition of the 3-Riemann tensor in 
terms of the 3-Christoffel symbols it is easy to obtain the 
linearized approximation for the 3-Ricci tensor and the 
3-curvature respectively: 



Ri 



e 3 R 



y 1 



Ri 



e 3 R 



f(j h i)i 



'Ah, 



5 i? 



3 R := dlh 



kl IL kl 



'Ah 



(A26) 

(A27) 

(A28) 
(A29) 



where 3 A := d\ stands for the Euclidean 3-Laplacian. 
Concerning the sources (matter fields) we assume 



E 
Ji 



Eo 



- eE 
el 



eSij 



(A30) 
(A31) 
(A32) 



The Hamiltonian constraint Eq. (A7) when linearized 
reads 



or explicitly 



5 i? = WttGqE , 



dlihi - 3 Ah = l6irG Q E 



(A33) 
(A34) 



T he linearized version of the momentum constraint Eq. 
(|A|), reads 



8iK l { - diK = 8nG Ji 



(A35) 



Finally the dynamic Einstein Eqs. (AE) when lin- 
earized read 



d t K\ + d l djN - 3 R* , = 4ttG 



(S - E)S i i - 2S\ 



The linear limit of evolution Eq. ( All ) reads 
d t K + 3 AN = AttGo \S + E 



(A36) 



(A37) 



The self-consistency of the 3+1 equations up to first 
order imply that the zero order source fields must vanish 
identically: 



E o = = ^ = S°j 



(A38) 



The above equations (A.34), (A35) and (A36) are the 
3+1 decomposition of the 4+0 equations fl28|). 
The 3+1 splitting of the perturbations (p6|) is 



(A39) 
(A40) 



-4k0 , (A42) 
The 3+1 splitting of the Lorentz gauge Eq. (|29| ) 
d»% v = 

d t h - 2 K d t 4> + 2d,N l - 2d t N = for v = 
d t N l + djhij - diN - \dih + nd^ = for v = i 

(A43) 







7oo = 




7oi = 


-Ni 


lij = 


hij — 


7 = 


-(s 



Using Eq. (A25) the zero component of Eq. (A42) reads 
d t N + nd t 4> = -K . (A44) 



From Eqs. (A44) and (A37) one obtains a wave equation 
for N: 



-df t N + 3 AA> = Kd^4> + 4ttG 



S + E 



(A45) 
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On the other hand, differe ntiati ng the spatial components 
of the Lorentz gauge Eq. ( A43 ) w ith re spect to time and 
using consecutively Eqs. (A2S), ( |A25| ) and the momen- 
tum constraint Eq. (A35) one obtains a wave equation 
for N i : 



-d^ t N l + 3 AN' 1 = 2nd 2 ti $ - IGttGoP 



(A46) 



Then Eqs. (|A45[), (|A46[), and (|A47|) read respectively 



d? t N + 3 AN = 4tt— ( £ matt + 5 mal 1 



F 

+ ^ + ^) d ^-§r 3 A4> . (A56) 



Finally, the linea rized dynamic Einstein Eqs. ( |A3q ) to- 
gether with Eq. (A22) and the spatial components of the 
Lorentz gauge Eq. ( A43 ) lead to a wave equation for hij 



-d} t K 3 + 3 Ah l3 = -2 K df 3 c 



inG 



(S - E)S i j - 2SV 
(A47) 



From Eq. (A47) one obtains 



d? t h + 3 Ah = -2k 3 A<j) + 8nG 



S-3E 



(A48) 



Another wave equation which is not independent from 
the above is obtained from the linearized Hamiltonian 



const raint Eq. (A34) and the Lorentz-gauge conditions 
Eq. ( |A43l) : 



-iff, -+ ( ^h-N) = -«(flg + a A).>- \{^G„E 



(A49) 



which is obtained by combining Eqs. ( A45 ) and (|A4 
This can be written as 



(-dl + 3 A) flh -K0-N 



2n 3 Acj>- 16ttGqE 



The combin ation of Eqs. flA45|) , ( |A47|) and (|A48j ) as 
given by Eq. (A41) provides the wave equation 



(-d 2 tt + 3 a) % = -iKdl 



XQitG^Sij . 
(A51) 



Therefore Eqs. ( |A46D , ( |A50| ) and ( |A5l| ) recover the 4+0 
wave Eq. (|3fj) in the Lorentz gauge. 

Now according to the 3+1 splitting of the energy- 
momentum tensor, and using Eq. (^Tj) we have 



-^+^ = -16^4^ 



(A57) 



d 2 h ■ 



3 Ah ll = 8tt^ 
F 



(q jTi \ri oomatt 

l^matt — -G'matt )0 j — ^O ij - 



(A58) 



From Eq. ( A48 ) one obtains 
Go 



d?~h + 3 Ah = 8tt 



F 



5E n 



2 ! K + f ) :a " 



3^D,0 
(A59) 



Moreover, from Eq. (A50) we get 



(-dl + 3 A) f \h - k4> - n) = -16n^E niat> 



Fn 



(A60) 



(A50) Whereas the Eq. (|A5l|) leads to 



(-dl + 3 A) % = -167r^£ matt 



F 



-2[k+^) {df^-5^ 



EL 

F 



(A61) 



Therefore Eqs. (A57), ( A60| ) , and (A61) are equivalent 
to Eq. ( |36| ) . In above equations one can employ the con- 
venient choice k = —Fq/Fq used in the 4+0 formulation 
to simplify the expressions. 



E — T n n — 



E n 



F 8irF G Q 



K 3 A0 , 



Tm — 



o^matt 
1 0i 

Fn 



(A52) 
(A53) 



Sij — Tij — 



T^matt 



S = T, = 



T 



matt I 



Fn 



8wF G t 

i r^(<9^-%M (A54) 
F' 

" 2 3 A<f> 



F ' 8^F G V « 



8nF, 



oG \ 6dt 



(A55) 



2. The weak- field approximation 

As in the 4+0 formulation, we consider slow vary- 
ing sources and neglect the time-derivatives, and also 
Smatt = P , S matt < p, J4 att « pc. With such con- 
siderations Eq. (A56) reads, 



'A N 



F' ~ 
— 



Air— p. 



(A62) 
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As in the 4+0 formulation 7^ = const. = N l are the 
regular solutions and the constants can be gauged out. 



44|) we can then obtain the so- 
Therefore the two non-trivial 



Combining ( |A62j ) and 
lution given by Eq. (p2 
degrees of freedom are N and tfi ( N l — 0, and, hij is 
given in terms of N and </>). The perturbed quantity N 
has the direct inte rpret ation of the Newtonian potential 
[cf. Eqs.([48|) and ( |A16| )]. In particular, for the massless 
case mo = 0, it yields, 



AN = 4tt 



GqP 
F 



F' 



ZGqFq 



(A63) 



where Eq. (j44|) was used. In the spherical symmetric 
case, ( |nj ) provides the exterior solution for N. 

Concerning the mass issue, let us focus on the lin- 
earized Hamiltonian constraint Eq. (A34). One can in- 
tegrate over a coordinate volume and use Gauss theorem 
to obtain 



16tt 



dihki - d k h)dS k = / Ed 3 x 



(A64) 



where the integrals extend to spatial infinite. The left- 
hand side corresponds precisely to the weak-field version 
of the ADM-mass formula [cf. Eq. (11.2.14) of Ref. @], 



M A 



DM 



1 

16tt 



(dih kl - d k h)ds k 



(A65) 



while the right-hand side of Eq. (A64) is the weak- field 
approximation of the Komar-mass formula. 



In the massless case mo = 0,from (A52) and (|44|), it 
turns 



E = 



1 - 



_2_(£o) 
2Go Fn 



(A66) 



Therefore, from Eqs. flAgg ) and ( |A66| ), one recovers 
Eq.(^p|) for the ADM-mass in the weak- field limit . 
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